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We review a series of works where the fundamental master equation is used to develop a micro- 
scopical description of evolution of non-equilibrium atomic distributions in alloys. We describe exact 
equations for temporal evolution of local concentrations and their correlators as well as approximate 
methods to treat these equations, such as the kinetic mean-field and the kinetic cluster methods. 
We also describe an application of these methods to studies of kinetics of Llo-type orderings in FCC 
alloys which reveal a number of peculiar microstructural effects, many of them agreeing well with 
experimental observations. 

I. INTRODUCTION 

Problems of evolution of non-equilibrium statistical systems attract attention in many areas of physics. These 
problems are of particular interest for configurational alloy kinetics — the evolution of the atomic distribution in 
non-equilibrium alloys. The microstructure and macroscopic properties of such alloys, e.g. strength and plasticity, 
depend crucially on their thermal and mechanical history — for example, on the kinetic path taken during phase 
transformations. Theoretical treatments of these problems usually employ either Monte G &f lo simulation — see e.g.EJ — 
or phenomenological kinetic equations for local concentrations and order parametersou. However, Monte Carlo 
studies in this field are time consuming, and until now they provided a limited information on the details of the 
microstructural evolution. Use of the phenomenological kinetic equations is more feasible, and Khachaturyan and co- 
workersoQ used this approach as a basis for discussing many microstructural effects. However, the phenomenological 
approach employs a number of unclear approximations — in particular, the extrapolation of linear Onsager equations 
for weakly nonequilibrium states to the nonlinear region of states far from equilibrium, and the relation between the 
phenomenological and microscopic approaches is often unclear. 

A consistent description of non-equilibrium alloys can be based on the fundamental master equation for the prob- 
abilities of various atomic distributions over lattice sites. The idea to employ this equation for studies of phase 
transformations was first suggested by MartinO. n-?°f ^ ne ^ as ^ decade this approach has been formulated in terms 
of both approximate and exact kinetic equationd9ii3 and was applied to many concrete problemsoEx In this pa- 
per we describe the main ideas and methods of this approach and illustrate them with an application to studies of 
microstructural evolution under Llo (CuAu I)-type orderings in FCC alloys. 

II. EXACT RELATIONS 

Following Ref. ^| we consider a substitutional alloy that includes atoms of m species p = Pi,P2> • • • Pm, in particular, 
vacancies for which p=v. Various distributions of atoms over lattice sites i are described by the different occupation 
number sets {n P i} where the operator n p i is unity when the site i is occupied by a p-species atom and zero otherwise. 
For each i the sum of operators n p i over all species p is unity, thus only m — 1 of them are independent, ft is 
convenient to mark the independent operators with special symbols, e.g. with greek letters: (n P i)i n d e p = n a i, while 
the rest operator denoted as Tiyi IS expressed via n a c 

n ri = \-y j n ai . (1) 

a 

The configurationally dependent part of energy H l can be written as 

Rt = E i r"p<"^ • E ^v"q^r fe + . . . (2) 

pq,i>j pqr,i>j>fc 

After elimination of the operators n r i according to ([!]) Eq. (Q) yields the interaction Hamiltonian H in terms of only 
independent operators n a i\ 
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H = ^ v°fn ai npj + + ^ v^riainpjn^k + ... (3) 

where effective interactions v" "f are linearly expressed via Vf"j q in (||). 

The fundamental master equation for the probability P to find the occupation number set {n a i} = £ is 

dP(0/dt = Y^[W(Z, V)P(V) - W(rj, 0P(0] (4) 
v 

where W(t;,r]) is the rj — ► £ transition probability per unit time. Adopting for this probabilitiy the conventional 
"thermally activated atomic exchange model"ul3, we can express W(£,rf) in (Q) in terms of the probabilitiy Wf q of 
an inter-site atomic exchange ("jump") qj <-> pi per unit time: 

W™ = n pl n qj Lo™ eM-PiKiw ~ 4?,<u )] = WWi? «p(^,«u ) ■ ( 5 ) 

Here uif q and E^ ^ are the "attempt frequency" and the "saddle point energy" assumed to be independent of alloy 
configuration; (3 — l/T is the reciprocal temperature; is the initial (before the jump) configurational energy of 

jumping atoms, and 7^ q = (x>^ q exp[— f3{E pi qj )] is the configurationally independent factor in the jump probability. 
If we accept for simplicity the pair interaction model, i.e. retain only the first term in Eq. (0), then the operator 
i?p™ qj - in (|J) may be expressed in terms of formal variational derivatives of the hamiltonian ([|) over n P i and 
H< = SHt/Snpi and ff* iqj . = S 2 H l / '5n pl S nqj : 



n a 



where the last term corresponds to the substraction of the "double-counted" interaction between the jumping atoms. 
The employed neglection of a possible configurational dependence of 7? q in (j^) is actually not essential, and one can 
also use any form of WFj q obeying the detailed balance principle. 

It has been shown irJa that in studies of practically interesting problems the true vacancy-mediated atomic exchange 
mechanism can usually be replaced by some equivalent direct exchange model. For example, instead of a real binary 
alloy ABv with vacancies and the vacancy-mediated atomic exchanges A<->v and B<->v we can consider a more simple 
model of a binary alloy AB with the direct A<->B exchange and only one independent variable n^i = n>i in Eqs. (^)-@ 
for each site i. Discussing below for simplicity only the binary alloy case we can seek the distribution function P(£) 
in (H|) in the form of a "generalized Gibbs distribution" : 

P{n i } = exp[/3(0 + ^A i n i -Q)]. (7) 

i 

Here the operator Q is an analogue of the Hamiltonian H in (|j|) : 

Q = ^ dijUiUj + ^jkriiUjUk + • • • ; (8) 

i>j i>j>k 

the "local chemical potentials" Xi and "quasi-interactions" ai...j (being, generally, both time and space dependent) 
are the parameters of the distribution; and the generalized grand canonical potential f2 = Q{Xi, g%...j} is determined 
by the normalizing condition: 

Q = -TlnTrcxp[/3(^ X.n, - Q)]. (9) 

i 

Multiplying equation (||) by operators rii, riirij, etc., and summing over all configurational states, i.e. over all number 
sets {rii}, we obtain the set of equations for the averages gij...k = (niUj . . . "Jr)i-p Tr(?iirtj . . . rikP), in particular, for 
the mean occupation c$ = (rii) = gi. After certain manipulations described inO'tj these equations can be written as: 

^ 9U...H = 7«« ( 2 s[nh ( D 7i) exp(£>+)n-nX- • ■ ■ n k) + {i —* j, . . . k}. (10) 

Here 7^ is 7^ B ; ri! i is 1 — nf, {i — ► j, . . . k} denotes the sum of expressions obtained from the first term by index 
permutation; the operators Df t are 
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Di = ^[(X s + (H-Q) s }±{s^i}; (11) 
and (H — Q)i is the variational derivative 5(H — Q)/Srii: 

(H - Q)i = { v ij-k ~ aij...k)nj ■ . -n k . (12) 

j>...>k 

Eqs. (Qjl^) enable us to derive the microscopic expression for the free energy F of the nonequilibrium stateBi: 

f = n + J2 x i c i+ < H -Q > ■ ( 13 ) 

i 

The function F — F{ci, 9i...j} obeys the generalized first law of thermodynamics, 

dF = Xjdcj + } {vi...j - a,i„.j)dgi,„j , (14) 

i i>...>j 

and has a fundamental property not to increase under spontaneous evolution, similarly to the Boltzmann's not 
decreasing entropy: 

dF/dt < 0. (15) 

The stationary state (being not necessarily uniform) corresponds to the minimum of F over its variables Ci and gi...j 
at the given N = J2i °i '■ 

dF/dci = \i = fi; dF/dgi...j = j - ai...j = (16) 



where fi is the Lagrange factor. Eqs. ( |16[ ) are the usual Gibbs relations for the parameters A.; and ai...j in the 
distribution <M) for the stationary case. 



III. KINETIC MEAN-FIELD AND KINETIC CLUSTER APPROXIMATIONS 

To approximately solve kinetic equations (|l^) one can use the regular approxi mate , methods of statistical physics, 
such as the mean- field approximation (MFA), the cluster variation method (CVM)EiK3, and also its simplified version, 
the cluster field method (CFM)EI. In both the kinetic MFA^nd kinetic CFM (KMFA and KCFM) the equations for 
Ci(t) are separated from those for gi...j(t) and have the forml3~1i3: 

dci/dt = 2 M H sinh[/?(Fj - Fj)/2] (17) 

3 

where F{ = Aj = dF/dci, while My = My{cfc} and F = F{ci} is the MFA or CFM expression for the generalised 
mobility and the free energy of a nonuniform alloy. For simplest approximations, such as the KMFA or the kinetic pair- 
cluster approximation (KPCA), these expressions can be written analytically. For example, the KMFA expressions 
for My and Fi in an alloy with only pair interactions V? q in (|^) are 

Fj = In — + y^fcCfc. (18) 

c * fe 

Here c- = 1 — c,, and u ik — V^ A — V® B is the "asymmetrical potential"!. Substituting; Eqs. ( |l8|) into ( |l7|) we obtain 
the analytical KMFA equation for Ci(t). A similar equation is obtained in the KPCAa. _ _ 

The usual phenomenological kinetic equations, in particular, those used by Khachaturyan and coworkersoH, corre- 
spond to the linearization of the KCFM equations ( |l7| ) in (Fi — Fj) and neglecting the c^-dependence in the mobility 
Mij. Such approach is usually sufficient for qualitative considerations, but for some-problems it can lead to a notable 
distortion of both the time scale and other details of the microstructural evolutionllJ. 

KMFA or KCPA are usually [Sufficient for studies of main kinetic features of spinodal decomposition, as well as 
orderings in the BCC latticeEalia. However, for more complex orderings, e.g. LI2 and Llo orderings in the FCC 
lattice, MFA and PCA are known to be insufficient and more precise methods are necessary, such as the CVM or 
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CFMcirES. Recently we suggested a simplified version of CVM, the tetrahedron cluster field method (TCFM)E3, that 
combines a high accuracy of CVM in describing thermodynamics with great simplification of calculations making it 
posible to develop its kinetic generalization, KTCFmH Similarly to the KMFA and KPCA, the KTCFM provides 
explicit equations for the mobility Mij and the local chemical potential Fi in Eq. (|l7j ) via mean occupations c&, and 
for each site i these equations can be reduced to a system of four nonlinear algebraic equations which can easily be 
solved using Newton's method. 

Let us now make remarks about effective interactions Vi...j in the Hamiltonian (j^) for real alloys. These interactions 
include the "chemical" contributions vf j which describe the energy changes under permutations of atoms A and B in 
the rigid lattice, and the "deformational" interactions vf • related to the lattice deformations under such permutations 
The chemical contributions are estimated from either first-principle calculations or fitting to some experimental dataE3, 
but for the long-ranged deformational interactions such methods can not be directly used. A microscopical model for v d 
in dilute alloys which includes only one experimental parameter was suggested by KhachaturyanH. The deformational 
interaction in concentrated alloys can lead to some new effects being absent in dilute alloys, for example, to the lattice 
symmetry changes under phase transformations, such as the tetragonal distortion under Llo ordering, and earlier 
these effects were described-only phenomenologicallyB. Recently we suggested-a microscopical model for calculations 
of v d in concentrated alloysO which generalizes the Khachaturyan's approach!. Unlike the case of dilute alloys, this 
deformational interaction turns out to be essentially non-parwise, and it includes two parameters which can be found 
from experimental data about the lattice distortion under phase transformations. 



IV. METHODS OF SIMULATION OF Ll ORDERING 



Up to recently most of theoretical treatments of kinetics of alloy ordering considered only simplest B2 (CuZn- 
type) orderings with just two types of antiphase-ordered domain (APD) and one type of antiphase boundary (APB) 
separating these APDs. Yet ordered structures in real alloys are usually more complex and include many.-types of 
APD. For example, under the DO3 (Fe3Al-type) ordering on the BCC lattice there are four types of APDEil, while 
under the Ll 2 (C^Au-type) or Ll ordering on the FCC lattice there are four or six types of APD, respectively. 
It results in a number of peculiar kinetic features that are absent for the simple B2 ordering. In Refs. |l8],[l0] we 
discussed such features for the DO3 and Ll2-type orderings. Below we consider the Llo-type orderings for which the 
microstructural evolution turns out to be still more complex and interesting. 

To study this evolution we made simulations of Al^Llo transformations after a quench of an alloy from the 
disordered Al phase to the single-phase Llo state. For these simulations we used the master equation approach 
and the KTCFM described above. We considered five alloy models with different types of chemical interaction: 
the second-neighbor (or "short-range") interaction models 1, 2, and 3 with v\ = 1000 K and the ratio ^2/^1 equal 
to (-0.125), (-0.25) and (-0.5l_respectively; the "intermediate-range" fourth-neighbor-interaction model 4 with v n 
estimated by Chassagne et alo from their experimental data for Ni-Al alloys: v\ = 1680 K, V2 = — 210 K, v% = 35 K, 
and i»4 = — 207 K; and the "extended-range" fourth-neighbor-interaction model 5 with v\ — 1000 K, V2/V1 — —0.5, 
V3/V1 = 0.25, and v^/vi = —0.125. The deformational interaction v d for all these models was estimated as described in 
section 3 with the use of experimental data for Co-Pt alloys. The critical temperature T c for our models corresponds to 
the stoichiometric composition c = 0.5 and is 614, 840, 1290, 1950 and 2280 K for model 1, 2, 3, 4 and 5, respectively. 

The distribution of mean occupations Cj under alloy ordering can be described in terms of both long-ranged and 
local order parameters. For the homogeneous LI2 or Llo ordering the distribution Cj = cCR.^) (where Ri is the FCC 
lattice vector) can be written in terms of three long-ranged order parameters f] a , see e. gJ3: 

c, = c + 771 exp(ikiRi) + r] 2 exp(ik 2 R,:) + rj 3 exp(ik 3 Ri) (19) 

where k Q is the superstructure vector corresponding to i] a : {ki, k2, k.3} — {[100], [010], [001]}2-7r/a. For the cubic LI2 
structure |^i| = lr/2 = [773!, 771772% > 0, and four types of ordered domain are possible. In the Llo-ordered structure 
with the tetragonal axis a, a single non-zero parameter rj a is present which is either positive or negative. Therefore, six 
types of ordered domains are possible with two types of APB. That separation of two APDs with the same tetragonal 
axis will be for brevity called the "shift- APB" , and that separation of the APDs with perpendicular tetragonal axes 
will be called the "flip- APB" . The transient partially ordered alloy states can be conveniently described in terms of 
local squared order parameters 77^ defined int3. The simulation results in figures 1-5 below are usually presented as 
the distributions of quantities rjj = 77^ + 77^ + 77^ (to be called the '^-representation" ) : the grey level linearly varies 
with ijf between its minimum and maximum values from completely dark to completely bright, and this distribution 
is similar to that observed in the transmission electron microscopy (TEM) imagesc3~L3. 

The simulations were performed in FCC simulation boxes of sizes VJ, = L 2 x H (where L and H are given in units 
of the lattice constant a) with periodic boundary conditions. We used both 3D simulations with H = L and quasi-2D 
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simulations with H = 1, and all significant features of evolution in both types of simulation were found to be similar. 
Each of figures 1-5 below includes all FCC lattice sites lying in two adjacent planes, z = and z — a/2, thus it shows 
4L 2 lattice sites. The initial as-quenched distribution Cj(0) was characterized by its mean value c and small random 
fluctuations 8ci] usually we used Sci — ±0.01. 

V. KINETICS OF A1->L1 TRANSFORMATION 

To avoid discussing the problems of nucleation, in this work we consider the transformation temperatures T lower 
than thfijOrdering spinodal temperature T s . Then the evolution under the Al— >Llo transition includes the following 
stagesErEj: 

(i) The initial stage of the formation of finest Llo-ordered domains when their tetragonal distortion still has little 
effect on the evolution and all six types of APD are present in microstructures in the same proportion. 

(ii) The imtermediate stage which corresponds to the so-called "tweed" contrast in TEM images. The tetragonal 
deformation of the Llo phase here leads to a predominance of the (llO)-oriented flip-APBs in the microstructures, 
but all six types of APD are still present in similar proportions. 

(iii) The final, "twin" stage when the tetragonal distortion of the Llo-ordered APDs becomes the main factor 
determining the evolution and leads to the formation of the (llO)-type oriented bands. Each band includes only two 
types of APD with the same tetragonal axis, and these axes in the adjacent bands are "twin" related, i.e. have the 
alternate (100) and (010) orientations for the given set of the (llO)-oriented bands. 

The thermodynamic driving force for the (100)-type orientation of flip-APBs is the gain in the elastic energy of 
the adiajeent APDs: at other orientations this energy increases under the growth of an APD proportionally to its 
volumealiil For an APD with the characteristic size I, surface Sd, tetragonal deformation e, and shear constant 
c s , this force begins to affect the microstructural evolution when the volume elastic energy E v el ~ c s e 2 lSd becomes 
comparable with the surface energy E s ~ aSd where a is the APB surface energy. The beginning of the tweed stage 
(ii) corresponds to the relation E v el ~ E s or to the characteristic APD size 

l Q ~ <?/c s e 2 . (20) 

The distortion e is proportional to the order parameter squarecQ and below it is characterized by its maximum value 
Era corresponding to T — and c = 0.5. 

Some results of our simulations are presented in figures |l|-^. The symbol A or A, B or B and C or C in these 
figures corresponds to an Llo-ordered domain with the tetragonal axis along (100), (010) and (001) and the positive 
or negative value of the order parameter 771 , 772 and 773 , respectively. Frame EH is shown in the /^-representation: the 
grey level linearly varies with 77^ between its minimumr-and, maximum from completely dark to completely bright, 
which corresponds to the usual bright-field TEM imagesc3ca 



5 



FIG. 1. Temporal evolution of the extended-range-interaction model 5 under the phase transformation Al— >Llo shown in 
the ^-representation for the simulation box size Vb = 128 2 x 1 at c — 0.5, the reduced temperature T 1 = T/T c — 0.7, e m | = 0.1 
and the following values of the reduced time t' = t<y i3 : (a) 10; (b) 20; (c) 50; (d) 100; (e) 250; and (f) 280. The symbol A, A, 
B, B, C or C indicates the type of the ordered domain as described in the text. The single, the double and the thick arrow 
indicates the splitting APB process, the quadruple junction of APDs, and the fusion-of-domain process, respectively, mentioned 
in the text. 
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Figures 1-5 illustrate quasi-2D simulations for which microstructures include only the edge-on APBs normal to the 
(001) plane. The above-mentioned elimination of the volume-dependent elastic energy in such geometry is possible 
only for the (100) and (OlO)-ordered APDs A or A and B or B separated by the (llO)-oriented flip-APB, while in 
the (OOl)-oriented domains C and C this elastic energy is always present. Therefore, the tweed stage (ii) in these 
simulations corresponds to both the predominance of (llO)-oriented APBs and the decrease of the number of domains 
C and C in the microstructures. 




FIG. 3. As figure g, but for the intermediate-range-interaction model 4 at T' = 0.67 and the following values of t': (a) 10; 
(b) 20; (c) 50; (d) 170; (e) 200; and (f) 700. 



Let us hrst consider figures |]-|| in which frames [l]a-[i]b, ||a and correspond to the initial stage; frames [T]c— |T|ci, 
|b-f|;, and |b-|c, to the tweed stage; and the restirames, to the twin stage. At both the initial and the tweed stage 
we can observe the following features of evolutior£9: 

(a) The presence of abundant processes of fusion of in-phase domains which are one of the main mechanisms of 
domain growth at these stages. 

(b) The presence of peculiar long-living configurations, the quadruple junctions of APDs (4-junctions) of the type 
A1A2A1A3 where A2 and A3 can correspond to any two of four types of APD different from Ai and Ai. 

(c) The presence of many processes of "splitting" of a shift- APB into two flip- APBs which is finished by either a 
fusion of in-phase domains mentioned in point (a) (s — > / process), or a formation of a 4-junction mentioned in point 
(b) (s — ► Aj process). 

For example, s — > / processes can be followed in frames |l|a-|l|b; [j]c-[l]d; |l|d-l]e; |^c-^d; etc. The-fusion with the 
disappearance of an intermediate APD which initially separates two in-phase domains to be fusedE3 can be seen in 
the lower right part of frames [l]a-|l]b. Several long-living 4-junctions are seen in frames [l]a-[l]d and ^c-||d; and an 
s — > Aj process can be followed in the lower right part of frames [l]a-[l]c. Let us also note that the microstructural 
features (b) and (c) can be naturally explained by a significant excess of the surface energy of shift- APBs with respect 
to flip- APBs found in our CFM calculations for the systems under consideration-.— 

Frames ||a-[3|c! also display some (lOO)-priented and thin "conservative" APBsEHfEj. Such APBs are most typical of 
the short-range-interaction systems — seeL9 and figure ^ below — where they have a low surface energy (being zero for 
the nearest-neighbor interaction model). Under an increase of the interaction range (as well as temperature or the 
deviation from stoichiometric composition 5c = 0.5 — c) the anisotropy of this surface energy decreases, and so for 
model 4 such APBs are few, while for the extended-range-interaction model 5 they are absent at all. 

Frames [^c-lljd^ ||b-|]c, and ||b-|]c (as well as f§a-f§b)_ illustrate a (HO)-type alignment of APBs between APDs A 
or A and B or B and a "dying out" of APDs C and C at the tweed stage. Frames |l|c and [j]d also show that in the 



7 



simulation with a realistic distortion parameter |e m | = 0.1 (fitted to the structural data for CoPt) the APD size Iq 
( f2p| ) characteristic of the tweed stega is about (20 — 40) a. It agrees with the order of magnitude of Iq observed in the 
CoPt-type alloys FePt and FePdNll 

Discussing the final, twin stage of the evolution we fiijsL-paention some characteristic configurations observed in 
experimental studies of transient twinned microstructuresclc3c3: 

(1) "semi-loop" shift-APBs adjacent to the twin band boundaries; 

(2) "S-shaped" shift-APBs stretching across the twin band; 

(3) short and narrow twin bands — usually with one or two shift-APBs near their edges-lying within the main twin 
bands; 

(4) an alignment of shift-APBs in the final, "nearly equilibrium" twin bands: within a (100) oriented band in a 
(llO)-type polytwin the APBs tend to align normally to sopfj-direction n = (cos a, sin a, 0) with a "tilting" angle a 
which is less than 7r/4 in CoPt and is close to zero in CuAuc3i23. 

Comparing our results with experiments one should consider that due to the limited size of the simulation box the 
twin band width d in our simulations has the same order of piaanitude as the APD size Iq ( p0| ) characteristic of the 
tweed stage, while in experiments d usually much exceeds Z(j£3~c3. Therefore, the distribution of shift-APBs within 
twin bands in our simulation is usually much more close to equilibrium than in experiments. In spite of this difference, 
the simulations reproduce not only "nearly equilibrium" configurations (4) but also transient configurations (l)-(3) 
and elucidate their formation mechanisms. In particular, both the "semi-loop" and "S-shaped" shift-APBs are formed 
from regular-shaped approximately quadrangular APDs (characteristic of early twin stages) due to the disappearance 
of adjacent APDs which are "wrongly-oriented" with respect to the given twin band; it is seen, for example, in frames 
|d-p, |e-|f, |c-|e, etc. The formation of short and narrow twin bands with the edges touched by shift-APBs is 
illustrated by frame |^d (which is strikingly similar to the experimental microstructure shown in Fig. 2 of Ref. |3|) and 
also by frame ||d. 

The alignment of shift-APBs mentioned in point (4) is illustrated by frames ^f and ^f, as well as |^d and ||a-^d. 
These frames show that the tilting angle a sharply depends on the interaction type, particularly on the interaction 
range, as well as on the concentration c and temperature T. In particular, for the extended-range-interaction model 
5 this angle is close to 7r/4; for the intermediate-range-interaction model 4 (which seems to more realistically describe 
properties of CoPt-type alloys) the angle a is less than tt/4, i.e. APB planes are tilted to the tetragonality axis; 
and for the short-range-interaction models 1 and 2 the APB planes tend to be parallel with the tetragonality axis, i. 
e. a ~ 0. Such interaction-dependent alignment of shift-APBs can be explainedEj by the competition between the 
anisotropy of their surface energy — which for both the intermediate and short-ranged-interaction systems tends to 
orient APBs parallel with the tetragonality axis decreasing the angle a, and a tendency to minimize the APB area 
within the given twin band which, corresponds to a — 7r/4. Therefore, the comparison of experimental tilting angles 
with the theoretical calculationscJ can provide both qualitative and quantitative information about the interatomic 
interactions in an alloy. 
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FIG. 4. As figure | 
60; and (d) 120. 



but for the short-range- interaction model 1 at T' = 0.9 and the following values of t'\ (a) 30; (b) 40; (c) 



ture |] illustrates evolution for model 1 which describes the short-ranged-interaction systems such as alloys Cu- 
. . The microstructures for such systems include many conservative APBs mentioned above, and the shift- APBs in 
the final frame |Jd are "step-like" consisting of (100)-type oriented conservative segments and small non-conservative 
ledges (being similar to the APBs observed in the L 12-ordered CU3A.U alloyed). These step-like APBs can be viewed as 
a "facetted" version of tilted APBs seen in frames |]f, |]f, ||c and ||d. As mentioned, under an increase of temperature 
T or "non-stoichiometry" 8c = 0.5 — c the anisotropy in the APB energy rapidly decreases. It results in sharp, phase- 
transition- like changes in morphology of aligned APBs, from the "faceting" to the "tilting" , which is illustrated by 
frames [|a-||d. These morphological changes are realized via some local bends of facetted APBs illustrated by frames 
|s|a |^b. Therefore, this "morphological phase transition" is actually smeared over some intervals of temperature or 
concentration, but frames |^a-|^d show that these intervals can be relatively narrow. 







FIG. 5. As figure | but for model 2 at [e m | = 0.1 and the following values of c, T and t': (a) c = 0.5, T' = 0.77, t' = 350; 
(b) c = 0.46, T' = 0.77, £' = 350; (c) c = 0.5, T' = 0.95, t' = 300; and (d) c = 0.44, T' = 0.77, t' = 300. 

Finally, let us make a general remark about the kinetics of multivariant orderings in alloys, such as the DO3, LI2 
and Llo orderings considered in Refs. [Ic|JT^ -po| and in this work. It is well known that the thermodynamic behavior 
of different systems under various phase transitions reveals features of universality and insensitivity to microscopical 
details of structure, particularly in the critical region near thermodynamic instability points. The results of this and 
other studies of multivariant oderings show that such universality does not seem to hold for their phase transformation 
kinetics, anyway outside the critical region (which for such orderings is usually either quite narrow or absent at all). 
The microstructural evolution reveals a great variety of peculiar features, the detailed form of which sharply depends 
on the type of interatomic interaction, structure, the degree of deviation from stoichiometric composition, and other 
"non- universal" characteristics. 
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